Quasiparticle density of states in Bi2Sr2CaCu208+5 single crystals probed using 
tunneling spectroscopy at ultra-low temperatures in high magnetic fields. 
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Break-junction tunneling spectroscopy at temperatures 30-50 mK in high magnetic field is used 
to directly probe the quasiparticle density of states within the energy gap in a single crystal Bi2212 
high- Tc superconductor. The measured tunneling conductances dI/dV{V) in the subgap region 
have a zero flat region with no evidence for a linear increase of the density of states with voltage. A 
number of tunnel break-junctions exhibited dl /dViV) curves with a second energy gap structure at 
the average magnitude 2Ap_p/e = 13 mV. Our data cannot be explained by either a pure s pairing 
or a pure d^2_y2 pairing. 

PACS numbers: 74.72.Hs, 74.50.-|-r 



I. INTRODUCTION 



A careful investigation of the density of states (DOS) 
for quasiparticle excitations in high-Tc copper oxide su- 
perconductors is essential for the understanding of the 
mechanism responsible for superconductivity. Tunneling 
spectroscopy has played an important role in verifying 
BCS theory in conventional superconductors. The con- 
ductance of a tunnel junction is directly proportional to 
the quasiparticle density of states. The fine structure 
found in the tunneling conductance at voltages above 
the energy gap is a direct proof that the electron-phonon 
interaction is the coupling mechanism for superconduc- 
tivity in conventional superconductors. It is widely ac- 
cepted that the high-Tc superconductors (HTS) have 
dx^_y2 symmetry of the order parameter. Among the 
HTS, BiaSraCaCuzOg-i-rf (Bi2212) has often been stud- 
ied by the tunneling method. Although the tunneling 
characteristics obtained with Bi2212 single crystals show 
features related to the superconducting energy gap, the 
shape of the tunneling conductance dl/dV vs V is far 
from the ideal BCS density of states. The conductance 
curves reveal a strong broadening of the superconducting- 
gap structure, with a nonzero contribution at zero voltage 
and a linear increase in the subgap region (cusp-like fea- 
ture) (see e.g. Ref lHH . In addition, an unusual feature 
frequently observed in the tunneling data of Bi2212 is the 
"dip-hump" structure beyond the gap edge which is close 
to the resonant spin excitation energy^ 

In superconducting tunnel junctions, the difference be- 
tween d- and s-wave symmetry leads to a significant 
change in the current-voltage {I{V)) characteristics. For 
an intrinsic unshunted s-wave superconductor-insulator- 
superconductor (SIS) junction, the zero-bias conductance 
should be exponentially small at low temperatures (espe- 
cially in case of superconductors with high Tc). In con- 
trast, if the order parameter has nodes and dx2_y2 sym- 
metry, the I(V) characteristics of the d-wave SIS junc- 
tion in the vicinity of the zero bias should have a finite 



slope and thus a nonzero conductance at zero voltage 
(see e.g. Ref. 0). Experimental tunneling characteris- 
tics of d-wave SIS junctions do not agree with the the- 
oretical predictions. Won and Makii found that at low 
temperatures the subgap tunnehng conductance dl/dV 
should increase almost quadratically with V. Whereas, 
Yamada and Suzuki^ have recently shown that for co- 
herent tunneling, the tunneling conductance of a d-wave 
SIS junction is distinctly different when compared to in- 
coherent tunneling. The subgap region of the coherent 
SIS tunneling conductance is almost linear in V. They 
find that the quasiparticle tunneling in a Bi2212 mesa is 
mostly coherent. In general, the symmetry and pairing 
mechanism of the superconducting state remains contro- 
versial. For this reason, it is important to measure the 
subgap tunneling DOS of SIS Bi2212 junctions at very 
low temperatures and at high magnetic fields. While the 
results of measurements of the superconducting energy 
gap in SIS junctions are insensitive to thermal broaden- 
ing, measurement in the mK region are required to rule 
out thermal excitations as the origin of both the excess 
conductance at zero voltage and the broadening of the 
conductance peak. 

Magnetic field is essential, firstly, in order to suppress 
the Josephson current, which makes the measurement 
of the subgap tunneling conductance difficult. Secondly, 
magnetic field is very useful to distinguish relevant infor- 
mation in tunneling data from anomalous features due 
to critical current effects in weak links of the junction 
and the Bi2212 single crystal. Finally, it is necessary to 
study the effect of the magnetic field on the subgap and 
gap structures in the tunneling spectra. 

Over the last few years the local DOS in Bi2212 has of- 
ten been studied spectroscopically using a scanning tun- 
neling microscope (STM). During the process of the me- 
chanical cleaning, the crystal breaks between the BiO 
planes of two adjacent half-unit cells.® In the case of tun- 
neling measurements performed on the BiO plane, the 
large separation between the tip and the pair of Cu02 
planes (greater than 10 A) precludes the possibility of 
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tunneling directly into these planes. In this situation, the 
electronic states of planes other than Cu02 play a role in 
the tunneling process and must be taken into accountii 
In contrast to STM, in a break-junction fabricated on 
Bi2212 single crystals by using a precision setup, the tun- 
neling can occur along Cu02 planes because many plane 
edges are created. It should be pointed out that since 
Bi2212 single crystals can be easily cleaved, the possi- 
bility exists that the crystal shears along a plane rather 
than breaking cleanly, leading to the formation of a c- 
axis junction, but we will argue later that the tunneling 
in our break junctions is most likely to be in the afc-plane. 
There are a number of articles that have shown that a me- 
chanically controllable break juncti on is one of the better 
tunnel systems (see e.g. references 3' 8*^. Breaking the 
crystal at cryogenic temperatures and in high vacuum 
(or in an inert cryogenic fluid) guarantees two atomi- 
cally clean surfaces, thereby minimizing surface contam- 
ination. The distance between these surfaces can be me- 
chanically controlled. Increasing the distance results in 
a smaller junction, which is eventually reduced to only a 
few atoms. It is possible to form a tunnel junction with 
a vacuum barrier between the two foremost atoms. ^ Al- 
though, in view of the technical complexity, this method 
has not been extensively applied, it has provided reliable 
results so far, especially in the case of Bi2212 single crys- 
tals. At the present time this is the only method which 
allows a tunneling investigations of Bi2212 single crystals 
at ultra-low temperatures. In this paper, we describe the 
first experimental study of the tunneling DOS in copper 
oxide HTS at temperatures 30-50 mK in high magnetic 
fields up to 26 T using high-quality break-junctions fab- 
ricated on Bi2212 single crystals. These investigations 
are an extension of our previous tunneling studies of the 
Bi-compound HTSiSiii 



II. EXPERIMENT 

The Bi2212 single crystals were grown by a KCl- 
solution-melt free growth method. ^^■-'^^ It is known, 
that overdoping or underdoping of Bi2212 can be 
achieved by changes in the oxygen content or by cation 
substitutions j^^'^^ The first method requires an anneal- 
ing in oxygen or argon at very high temperatures. A 
careful characterization of the annealed samples reveals 
that changes in Tc are always accompanied by a severe 
degradation of the sample quality. Kinoda et ali^ have 
showed also that the annealing in oxygen can substan- 
tially increase the gap inhomogeneity. For this reason 
we have used only high quality as grown single crys- 
tals in which, the substitution of trivalent Bi for diva- 
lent Sr during growth reduces the hole concentration in 
the Cu02 planesii^ As the Bi/Sr ratio increases, the 
number of holes doped into the system decreases, which 
therefore pushes the system towards the hole-underdoped 
regime. Very recently Eisaki et ali^ have shown that as 
the Bi/Sr ratio in Bi2+a;Sr2_a:CaCu208+5 single crystals 




FIG. 1: The scanning electron micrograph of a crystal frag- 
ment where the composition measurement points are denoted 
by crosses 



approached 1, Tc increases from 82.4 K for a; = 0.2 to 
91.4 K for X = 0.07, 92.6 K for x = 0.04, and eventually 
to 94.0 K. Bi2212 single crystals with the excess Bi have 
repeatedly been used in tunneling study of temperature 
and doping dependence of the superconducting gap and 
pseudogap (see e.g. Ref. ITsI) . 

The quality of the crystals was verified by measure- 
ments of the dc resistance, ac susceptibility. X-ray diffrac- 
tion and scanning electron microscopy. The single crys- 
tal showed X-ray rocking curves with a width of about 
0.1° demonstrating the high quality and high homogene- 
ity of the samples. The values of the residual resistivity 
in our crystals were comparable with the lowest value 
previously reported for Bi2212. (see e.g. Ref. Il9|) . The 
composition of the crystals was studied using a Philips 
CM-30 electron microscopy with a Link analytical AN- 
95S energy dispersion X-ray spectrometer—. The actual 
cationic compositions of each crystal investigated were 
measured at several different places on the crystal and 
the scatter in the data was less than 2%. In Fig. ^ we 
show a scanning electron micrograph of a crystal frag- 
ment where the composition measurement points are de- 
noted by crosses. Previously, we have measured the Hall 
coefficient in several crystals and found the nearly linear 
relation between the excess Bi and the hole concentra- 
tion p. Our samples showed the well-known parabolic 
behavior for Tdp). 

Since the deviation from conventional superconductiv- 
ity should be most pronounced in the underdoped regime, 
for the present study, we have chosen three underdoped 
as grown single crystals Bi2.22Sri.55Cai.i7Cu2.oi08-i-5 
with Tc = 84 K and transition width ATc = 1.5 K. 



3 




50 100 150 200 250 300 
Temperature (K) 



FIG. 2: The temperature dependence of the in-plane re- 
sistivity pab for one of the investigating underdoped crystals 
(#1, p — 0.125). For comparison, we show also pab{T) for 
an overdoped {p = 0.2) single crystal with the same Tc = 84 
K (midpoint) as well as the data for a slightly underdoped 
sample with Tc = 87 K (p = 0.13) 



The same hole concentration for the investigated sam- 
ples p = 0.125 has been obtained from the Bi/Sr ra- 
tio as well as from the empirical relation Tc/Tc^max = 
1 — 82.6(p — 0.16)^ which is satisfied for a number of the 
HTS.^^ We used the fact that the substitution of triva- 
lent Bi for divalent Sr in the Bi-compounds reduces the 
hole concentration in the Cu02 planesii^ 

In order to demonstrate that the investigated samples 
are underdoped, in Fig. [21 we show the temperature de- 
pendence of the in-plane resistivity pab for one of the 
investigating crystals (#1, p — 0.125). For comparison, 
we show also Pab{T) for the overdoped (p — 0.2) single 
crystal with the same Tc = 84 K (midpoint) as well as the 
data for the slightly underdoped sample with Tc = 87 K 
{p = 0.13). One can see that as in all Bi-compound HTS, 
the magnitude of Pab{T) increases with decreasing carrier 
concentration. A typical T-linear behavior and a slightly 
upward curvature of pab are also seen in the slightly un- 
derdoped and the overdoped samples. Whereas for the 
underdoped sample {p = 0.125), pab deviates from high- 
temperature T-linear behavior at a characteristic temper- 
ature T* (indicated by the arrow in Fig . , as would be 
expected. This temperature in Ref. [2^ was identified as 
the pseudogap closing (opening) temperature T*. This is 
a further proof that we are measuring truly underdoped 
samples. 

The dimensions of the investigated crystals were ~ 
1 mm X (0.5 — 1) mm x (1 — 3) /xm. A four-probe 
contact configuration, with symmetrical positions of the 
low-resistance contacts (< 1J7) on both a6-surfaces of 
the sample was used. For the tunneling measurements, 
the sample was cooled to 30 mK at zero magnetic field. 
The tunnel junction was then fabricated in situ using a 
mechanical break junction technique. With the geom- 



etry of setup used, the break junctions forms so that 
tunneling occurs in the a6-plane. The I{V) characteris- 
tics and differential conductances dl/dV as function of 
V were measured using phase-sensitive detection tech- 
niques. We were able to fabricate a large number of tun- 
nel break-junctions at different places along the initial 
break with resistances from 100 to 120 kfi (for bias 
voltages around 300 mV) . Each junction was stable with 
reproducible characteristics at different magnetic fields. 
With our break-junction setup the magnetic field was 
oriented parallel to the Cu02 planes to an estimated ac- 
curacy of W l°|i& 

In a previous investigation we have measured break 
junction tunneling together with the a&-plane and out- 
of-plane resistivities for the same Bi2212 crystal at dif- 
ferent temperatures (4.2 - 250 K) for magnetic fields up 
to 20 T oriented parallel to the c-axis.^° In a magnetic 
field pab{T) reveals a significant increase in the width 
of the superconducting transition, while the value Tc for 
the onset of superconductivity remained constant. The 
data for PciT) shows evidence for a strong suppression 
of the superconductivity along the c-axis in a magnetic 
field of 20 T. The observed decrease of the Tc value by 
a factor 2 indicates a significant reduction of the energy 
gap. In contrast, for all tunnel junctions investigated 
the energy gap is almost independent of magnetic field 
(i/||c-axis)< 20 T)}^ For this reason, we can formally 
state that the observed gap spectrum for our break junc- 
tions is related to the a6-plane superconducting energy 
gap and therefore to tunneling in the a5-plane. At the 
same time, we are unable to say anything about the tun- 
neling direction in the a6-plane itself. For example, it was 
shown in Ref. 23] that the shape of spectra measured us- 
ing point contact tunneling spectroscopy of YBa2Cu307 
within the a6-plane in the [100] and [110] directions were 
quite different although the extracted magnitudes of the 
energy gap were in close agreement. The exact knowledge 
of the tunneling direction is essential for point contact 
tunneling spectroscopy where the main contribution to 
the quasiparticle current is given by Andreev reflection, 
which strongly depends on the wave vector. However, 
this is not the case for tunnel junctions. Hoogenboom 
et al.^ have shown on the basis of a detailed analysis of 
a large number of tunneling and ARPES spectra that 
the tunneling spectroscopy probes states along the entire 
Fermi surface. 



III. BREAK JUNCTION TUNNELING 
SPECTROSCOPY 

FigEIa-c) shows representative I{V) curves together 
with the tunneling conductances data dI/dV{V) for a 
100 il-, 115 kn,-, and a 6 kfi tunnel junction, respec- 
tively, at T ~ 40 mK, for different magnetic fields (sam- 
ple #1). All the I{V) characteristics exhibit characteris- 
tic features typical for superconducting tunnel junctions, 
with a fiat region around zero bias, consistent with the 
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FIG. 3: I{V) characteristics for the (a) 100 fl-, (b) 115 kfi- 
, and (c) 6 kSl-resistance tunnel junctions at T = 30 and 
50 mK, in different magnetic fields (sample #1). The insets 
in (a) and the upper insets in (b,c) show the corresponding 
tunneling conductances data dI/dV{V) in the subgap region. 
The lower insets in (b,c) show the same I{V) characteristics 
in the gap region, (d-f) the data points (squares) represent 
the density of states N{E) obtained from the deconvolution 
of Eq.Q (see text) 



expected zero tunneling conductance (the insets in FigO 
(a-c)), and a well-defined sharp increase in the tunnel 
current near ±(45 — 65) mV connected with the super- 
conducting gap. At large bias voltages up to ±0.6 V, 
the I{V) characteristics are linear with small deviation 
(increasing conductance) as expected for junctions with 
a good tunnel barrier without current leakage.^** The 
I{V) characteristics for the low-resistance tunnel junc- 
tion in the subgap region (FigUJa)) are also typical for 
low-resistance Josephson junctions (the arrows indicate 
the direction of the bias current sweep during the mea- 
surements). The I{V) curves of this type were hysteretic 
(versus bias sweep direction) and have the slope of this 
switching. This is evident even from (Fig|3Ia)). In ac- 
cordance with the direction of the bias current sweep, 
the Josephson current in zero magnetic field equals -76 
and +120 /iA. In our Bi2212 break-junctions we have 
observed the expected Fraunhofcr-like dependence of the 
Josephson critical current as a function of magnetic field. 
The large period of Fraunhofer pattern points to the ex- 
cessively small size of the tunnel junctions. After the 
suppression of the Josephson current in a magnetic field, 
the tunneling conductances dI/dV{V) of these junctions 
have a well-defined large zero fiat region around zero bias 



voltage (the lower inset in Figl^Ja)). The high resistance 
junctions do not show a zero-bias peak even at _ff = T. 

It should be pointed out that Joule heating effects 
could be important at the ultra low temperatures, espe- 
cially for low impedance junctions. Previously we have 
studied Joule heating effects in break-junctions fabri- 
cated on Bi22iaia and BizSrzCuOe-i-s (312201)^5 single 
crystals. It was found that the shape of the tunneling 
conductance dI/dV{V) and the value of the zero-bias 
dl/dV are very sensitive to a heating of the crystal re- 
gion near the tunnel barrier. The dI/dV{V) characteris- 
tics in our data (FigUfa-c)) show no evidence whatsoever 
of heating. The linearity of the current-voltage charac- 
teristics at high bias voltages (up to 0.6 V) indicates that 
the heating caused by the current injection in our mea- 
surements is negligible. The shape of the tunneling con- 
ductance dI/dV{V) is unchanged despite a variation in 
the junction resistance of three orders of magnitude. 

The measured tunneling conductances in the subgap 
region at > |10| mV can be fitted quite well to curves 
which are quadratic in V. The tunneling conductances in 
case of the high-resistance tunnel junctions in Fig|3Ib,c) 
in the subgap region do not show a large zero fiat re- 
gion. We have similar characteristics to those shown 
in Fig|3Ia-c) for several tens of tunnel break-junctions 
formed on each of the three single crystals. We have 
never observed, either a linear increase of the tunnel- 
ing conductances with voltage in the subgap region, or 
a nonzero zero-bias conductance for temperatures ~ 50 
mK. The conductance spectra with a flat bottom seems 
at first sight to be consistent with s-wave rather than 
(i-wave superconductivity. It should be pointed that the 
nonlinear subgap tunneling conductance in the vicinity of 
the zero bias has been observed previously in Bi22132ii£ 
and Hg-based cuprates^SiSi at 4.2 K, but the magnitudes 
of the background conductances at zero voltage were siz- 
able, suggesting a significant density of states of quasi- 
particles at zero energy. 

Our results are nevertheless puzzling since the tunnel- 
ing conductances in the subgap region in s-wave SIS junc- 
tions at mK-temperatures should have the zero flat region 
practically up to the gap voltages ±2A/e, and exhibit a 
very sharp conductance peak structure at = ±2A/e. 
This is clearly not in agreement with our results. De- 
spite the fact that the measurements were carried out 
at very low temperatures, the gap structure is strongly 
smeared with a broad conductance peak, and the shape 
of the tunneling spectra are not consistent with the ex- 
pected tunneling conductance for s-wave superconduc- 
tors at temperatures near K. On the other hand, the 
quadratic in V behavior of the subgap tunneling conduc- 
tances observed in this work is in agreement with the 
conductance of the SIS junction calculated by Won and 
Maki* with d-wave order parameter, although the shape 
of our experimental curves at voltages above the conduc- 
tance peaks differs strongly from the theoretical curves. 
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IV. QUASIPARTICLE DENSITY OF STATES 

In order to extract the quasiparticle DOS, N{E), 
from the tunnehng conductances in Fig|3{a-c), we de- 
convoluted the expression for the normahzed tunnehng 
conductanceSi 



d 



dI/dV{V)s _ 

dI/dV{V)N deV 



N{E)N{E - eV)dE. 



(1) 



This expression suggests that a tunnehng matrix element 
which weights a particular momentum direction is a con- 
stant and the tunneling spectrum directly relates to the 
DOS in the Cu02 plane. Although in some articles the 
tunneling spectrum was considered to be suggestive of 
an anisotropic matrix element, in Ref. it was shown 
on the basis of a detailed analysis of the tunneling and 
ARPES spectra that the tunneling matrix element does 
not have a strong dependence on the wave vector. Tun- 
neling spectroscopy therefore equally probes states over 
the entire Fermi surface »^ 

An iterative procedure for the deconvolution together 
with the method used to obtain the normalized tunnel- 
ing conductance is described in Ref. 0. In Fig.|3fd-f) 
the data points (squares) represent the result of this de- 
convolution for three different tunneling junctions from 
FigEIa-c), respectively. The solid lines show the exper- 
imental conductance curves [dl / dV) s / {dl / dV) m com- 
pared with the calculated curves (circles) using the N{E) 
obtained from the deconvolution of Eq.QJ. The good 
agreement with the data in Fig.|3Jd-f) confirms the valid- 
ity of the iterative procedure used to obtain N{E). The 
shape of the N(E) showed in Figs|2Id-f) is fully compati- 
ble with a smeared BCS DOS for Bi2212i2, The obtained 
N{E) over the entire subgap region can be fitted quite 
well to curves which are quadratic in V (dashed lines) . A 
quadratic rather than linear increase of N(E) in the sub- 
gap region seems to be strong evidence against d-wave 
symmetry. 

Note that a quadratic or nonlinear in V subgap tunnel- 
ing conductance in Bi2212 has been observed previously 
in STM experiments, for example in Ref. |7l l2al29l |. Franz 
and Millis^ have shown that a reasonable fit to the data 
in Ref. assuming a DOS with a d-wave order pa- 

rameter could be made provided the tunneling matrix 
element is anisotropic. Assuming an isotropic (constant) 
matrix element leads to sharper V-shaped gap structure, 
inconsistent with the experimental data^ As previously 
noted, the tunneling matrix element is most likely in- 
dependent of the wave vector;^ which implies that the 
STM spectra give the quasiparticle DOS directly, and 
are therefore inconsistent with a pure c?-wave order pa- 
rameter, in agreement with our results. For this reason, 
we are forced to conclude that it is not possible to ex- 
plain our data in terms of either pure s pairing or pure 
dx2-y2 pairing. Our results lead to the conclusion that 
Bi2212, in all likelihood, has a, d + s mixed pair state, as 
has been repeatedly suggested in the literature (see e.g. 
RefUm. 
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FIG. 4: The effect of the magnetic field on the gap-like struc- 
ture (marked by arrows) in the tunneling conductance dl /dV 
vs. V within the superconducting gap for different break- 
junctions on three single crystals at T = 50 mK and T = 30 
mK. (a,b) - two tunnel break-junctions with resistance nearly 
1 kS7 fabricated with the same single crystal #7. The inset 
shows the full gap spectra at and 26 T. (c,d) - the subgap 
regions of the spectra on an enlarged scale for two single crys- 
tals #5 and ^1 (the resistances nearly 300 Q). For clarity, 
the curves in (b-d) have been shifted vertically with respect 
to the lower curves. 



V. SUB-GAP STRUCTURE 

In addition, a large number of the tunnel break- 
junctions exhibited dl /dV{V) curves with gap-like fea- 
tures within the gap at voltages ~ 6 — 20 mV. With 
increasing magnetic field this structure broadens, dimin- 
ishes in amplitude and shifts to lower voltages. In Fig01 
we show the subgap regions of the differential conduc- 
tances dl/dV as a function of the bias voltage V in dif- 
ferent magnetic fields at T ~ 50 mK. (a,b) - two tunnel 
break-junctions with resistance of nearly 1 kQ fabricated 
with the same single crystal #7 (the low- voltage gap-like 
structure is indicated by arrows). The inset in (a) shows 
the full gap spectra at and 26 T in order to compare 
the amplitude of the sub-gap features with the gap peak. 
FigEIc,d) shows the sub-gap regions of the spectra on an 
enlarged scale for two single crystals #5 and #1 (with 
resistances of around 300 fl). For clarity, the curves in 
FigEIb-d) have been shifted vertically with respect to 
the lower curves. 

It might be thought that the observed sub-gap struc- 
ture are low-energy features related to the zero-bias peak 
which is caused by cZ-wave Andreev bound states^^'^"* or 
the Josephson effect. However, in case of the higher- 
resistance break-junction (Fig^Ja)) the zero-bias peak is 
lacking at = T. Nevertheless, both the sub-gap struc- 
ture and its field dependence are seen. In Fig^b-d), the 
sub-gap gap-like structure is also clearly visible despite 
the fact that the zero-bias peak is suppressed completely 
by the magnetic field. This data show that the zero-bias 
peak and the sub-gap gap-like structure in our case have 
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distinctly different origins. 

Taking into account the magnetic field dependence of 
the gap-like structure observed here, we identify this fea- 
ture with a small energy gap, the average magnitude of 
which is 2Ap_p/e = 13 mV. Electronic structure cal- 
culations for Bi2212 have long predicted the presence 
of BiO-derived bands and a splitting of Cu02 bands at 
the Fermi energy. If we assume that the break-junction 
tunneling averages the DOS equally around the Fermi 
surfacCf^ then the observed tunneling conductance must 
involve two separate order parameters derived from dif- 
ferent bands. While the sub-gap structure might be con- 
nected with a proximity induced energy gap in the BiO 
layersj^ the possible existence of a second small gap 
(2A < 8 - 20 meV) in Bi2212 has been previously sug- 
gested by Shen et al^, Mahan^'^ and Mallet et alM in 
order to explain the gap anisotropy observed in ARPES 
and vacuum tunneling spectroscopy of Bi2212. It should 
be pointed out that the weak sub-gap structure in Bi2212 
has been occasionally seen in the vicinity of the zero bias 
in STM spectra at 4.2 KSiSa and break-junctions at 10 
Km^ It is also possible that the sub-gap structure de- 
tected by us at mK temperatures can manifest itself as 
a "linear" increase in the subgap region in the vicinity 
of the zero bias as often observed in the STM spectra of 
Bi2212 measured at T = 4.2 K. 



VI. SUPERCONDUCTING ENERGY GAP 

We now turn our attention to the magnitude of the 
superconducting energy gap in the studied crystals. For 
the SIS junctions studied here, the peak-to-peak sepa- 
ration of the two main maxima of the dl /dV{V) curves 
must correspond to 4Ap_p/e. As can be seen in Fig|3fa- 
c), this value is different for the different break-junctions 
formed using the same single crystal (#1) and equals 96, 
93, and 124 mV. The applied magnetic field H \\ ab- 
plane should decrease only slightly the peak-to-peak sep- 
aration in dI/dV{V) curves due to the very high paral- 
lel critical magnetic field in Bi2212. For sample #1, we 
find a surprisingly wide scatter of the size of the energy 
gap 2Ap_p/e ranging from 46 to 78 mV for the more 
than 50 different break-junctions investigated. Similarly, 
for two other single crystals (#5 and #7), we obtained 
2Ap_p/e = 50 - 65 mV and 2Ap„p/e = 48 - 72 mV, 
respectively. 

Miyakawa et alm^^ have found a strong monotonic de- 
pendence the superconducting energy gap of the Bi2212 
on the doping concentration. However, we can for- 
mally exclude macroscopic sample inhomogeneity, as the 
origin of the observed variation, since 2Ap_p/e= 46- 
78 meV would imply''" a spread of hole concentration 
p = 0.2 — 0.15. Using the empirical relation for Tc vs p, 
this would imply values of Tc ranging from 80 to 95 K. In 
addition, the chemical composition of the crystal, mea- 
sured at several different points, is found to be identical 
(FigQ. The crystal is of a very high quality judging from 
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FIG. 5: Differential conductances dl /dV as a function of the 
bias voltage V for tunnel break-junctions fabricated on single 
crystals #7 and #1 with double peaks in different magnetic 
fields at T ~ 50 mK. For clarity, the curves have been shifted 
vertically with respect to the lower curves 



the sharp superconducting transition width ATc = 1.5 
K (in the magnetization measurements), and the small 
rocking curve width 50 ~ 0.1°. Following reference KEi 
this strongly suggests that such a wide range of the gap 
values originates from the presence of different supercon- 
ducting regions on the scale of a few coherence lengths. 
The effects of a nanoscale chemical inhomogeneity on Tc 
of Bi2212 was recently studied in Ref. If the sample 
has a nanoscale disorder with multiple energy gaps, the 
convolution of several DOS gives a strongly smeared sin- 
gle conductance peak (Fig|2Ia-c)), nevertheless, the DOS 
after the deconvolution of the conductance can show an 
additional maximum in the peak region (Fig|3ff)). In in- 
dividual cases the break-junctions can show double peaks 
in the conductance at voltages nearly 2Ap_p/e (FigE)), 
that is indicative of two nanoscale superconducting re- 
gions in the region of the junction. 



VII. CONCLUSION 

In conclusion, we have studied the tunneling density 
of states using high-quality break-junctions fabricated 
on Bi2212 single crystals. Our experiments provide a 
crucial test for the presence of excitations within the 
superconducting energy gap. The measured tunneling 



7 



conductances in the subgap region show a zero flat re- 
gion. We have never observed either a hnear increase of 
the tunneUng conductances with voltage in a subgap re- 
gion or a nonzero zero-bias conductance at the temper- 
atures 30-50 mK. The deconvoluted DOS in the whole 
subgap region can be fitted quite well to curves which 
are quadratic in V. A quadratic rather than linear in- 
crease of N{E) in the subgap region seems to be strong 
evidence against pure dx2_y2 symmetry. The major part 
of the tunnel break-junctions exhibited dI/dV{V) curves 
with the second energy gap structure with an average 
magnitude 2Ap_p/e = 13 mV. We cannot explain our 



data with either pure s pairing or pure d pairing. The 
wide range of the superconducting gap values observed 
suggests that the results of investigations of the doping 
dependence of the superconducting gap in HTS by tun- 
neling spectroscopy need to be treated with care. 
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